Frequency selection by soliton excitation in nondegenerate intracavity downconversion 
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We show that soliton excitation in intracavity downconversion naturally selects a strictly defined 
frequency difference between the signal and idler fields. In particular, this phenomenon implies 
that if the signal has smaller losses than the idler then its frequency is pulled away from the cavity 
resonance and the idler frequency is pulled towards the resonance and vice versa. The frequency 
selection is shown to be closely linked with the relative energy balance between the idler and signal 
fields. 



Exchange of ideas between nonlinear optics and non- 
linear dynamics of spatially distributed systems has pro- 
duced a series of interesting results over the last decade, 
and has opened up one of the most active areas of current 
research. In this particular work we will consider novel 
phenomena associated with soliton excitation in the prac- 
tically and fundamentally important area of nondegener- 
ate intracavity down-conversion. 

The essence of parametric down-conversion is virtual 
absorption of a pump photon at frequency lu p with a sub- 
sequent reemission of two photons with frequencies 
and uj s , where indices i and s stand, respectively, for the 
idler and signal fields. Down-conversion can be realized 
both in free propagation and in intracavity schemes. The 
latter takes advantage of positive feedback provided by 
the mirrors and thus transforms the passive free propaga- 
tion scheme into an active generator or optical parametric 
oscillator (OPO) @,|. 

Influence of the transverse degrees of freedom on the 
quantum || and classical [p|-pd| aspects of the parametric 
processes have recently become a subject of significant ac- 
tivity. Among the main attraction points on the classical 
side are localised structures |)| [llj . Bright diffractionless 
localised excitations inside an optical cavity supported 
by different nonlinearities (cavity solitons) have been re- 
cently observed experimentally []12| and suggested for 
all-optical processing and storage of information |p|,|l3|, 
see also recent reviews JIJ]. The large quadratic nonlin- 
earities of artificially phase matched materials make 
parametric cavity solitons |s|-[ll| particularly attractive 
for practical application, especially where fast material 
response is an issue. 

Down-conversion processes can be divided into degen- 
erate and nondegenerate. In the former case idler and 
signal photons are identical while in the latter they differ 
in frequencies and/or polarizations. It has been shown 
that transverse patterns in nondegenerate OPO and 
soliton dynamics in nondegenerate free propagation pd| ] 
have qualitative differences from their degenerate coun- 
terparts. The difference can be formally identified as due 
to an additional symmetry, in the differential phase of 
signal and idler fields pd| . This symmetry is suppressed 
in the degenerate case. As a result the frequency of the 



signal component of any solution, including solitons, is 
exactly half the pump frequency, lu s = uji = lu p /2. On 
the other hand, in the nondegenerate case any arbitrary 
frequency difference 2f2 between the idler and signal fields 
still satisfies the condition oj p = uj s + LUi . This raises the 
question of whether there are any physical constraints on 

ft. 

While in free propagation the value of f2 is limited 
only by phase matching conditions, this problem becomes 
more subtle in the OPO because cavity effects come into 
play. A review of early works on this issue, exploring 
approaches based on the plane wave approximation, can 
be found in . More recently Longhi M has shown that 
if diffraction is included then fl becomes a function of 
the magnitude of the transverse component of the signal 
and idler wave vectors, k\ = — kj_, of the exact trav- 
elling wave solution which exists in this system Q. For 
the fixed OPO parameters {k^l can take any values from 
certain continuous bands and hence so can Q. 

The primary object of this Letter is to demonstrate 
that Q is constrained when a cavity soliton is excited in 
the nondegenerate OPO. We show how this follows from 
a general relation between the soliton energies and the 
cavity losses. Existence of this relation, which has not 
been previously identified in this context, seems to be 
closely related to survival of differential phase symmetry 
in the presence of cavity losses. We use this symmetry 
to derive aproximate formulae for fl in certain limits. 
Understanding this problem is not only practically im- 
portant, but also holds the key to construction of entire 
families of cavity solitons. 

Mean-field equations describing interaction of the sig- 
nal, idler and pump waves in OPO f|,|],|| can be pre- 
sented in the dimcnsionlcss quasi-Hamiltonian form 
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where H is the following functional: 

/ dx[-a s \d x E s \ 2 - ai\d x Ei\ 2 - a p \d x E p ^ 2 



t is 



H 

Ss\E s \ 2 

S^ 2 + S p \E p \ 2 + (E p E*E* + ^E* S E* + c.c„ 
the time measured in the units of r c 6, where r c is 
the cavity round-trip time and b is an arbitrary scal- 
ing constant. 7,„ = T m b/2, where T m are the effec- 
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tive mirror transmitivities. x — X[2fc s /(6L)] 1 / 2 and 
a m = k s /k mi where X is the transverse coordinate in 
physical units, k m are the longitudinal components of 
the wave vectors and L is a cavity round-trip length. 
5 m = bT c (uj m — w™ 11 ) are the detunings from the cavity 
resonances Physically, validity of these equations 

requires small losses and detunings: j m ,S m << bn. As- 
suming, for simplicity, that lu s is close to u>i, we can 
fix a St i — 2a p — 1, but still allow differences in 8 S , Si 
and 7 S , 7^. Then the phase matching conditions imply 
n(u> s ) — n(u>i) ~ n(bj p ), where n(u>) is the linear re- 
fractive index and the following scalings can be derived 
for the field variables, fi is proportional to the external 

pump field E p : E p = ^^2(6 2 + j 2 )/f UJp /(bT c xuJs), here 
X is the effective quadratic susceptibility. The physi- 
cal fields E m are given by = \/2E s ^e 1 ^^ 2 /(bT c x^ s ), 
£ p = 2e^(E p + nJsl+ll)/{bTcXUs), where = 

—atan(j p /Sp). 

We seek localized solutions of Eqs. (1) in the form 
E m (x,t) = A m (x, A)e in "*, where tl Sii = ±A ± (S s - 
5i)/2, flp = 0. Then A m obey the set of differential 
equations 

-ij s A s = (d% + 5- A)A S + (A p + fi)A*, 

-ijiAi = {d 2 + 5 + A)A l + {A p + M ) A* s , (2) 

-i2 lp A p = (d 2 + 25p)Ap + 2A s A l , 

where <5 = (5 S + 5i)/2. The previously introduced 
frequency difference 2f2 between the idler and signal 
fields is linked with the parameter A by the formula: 
2Q = (2A + S s — Si)/(bT c ) We are interested in bright 
single-hump solitons, implying A m (x — +00) = and 
d x A m (x — 0) = 0. Existence of such solitons in the 
parameter region where the trivial zero solution bifur- 
cates subcritically can be predicted for A = 0, j s = 7« 
by analogy with the well studied degenerate case, where 
this condition reads SS p > j s j p |]6|,p|Jl0| . In the nondegen- 
erate case solitons have been found as a result of direct 
numerical simulations of Eqs. (^) starting from either 
'random' [|| or localized Q initial conditions. Subcrit- 
ical bifurcation and the related phenomenon of optical 
bistability in the nondegenerate OPO have been demon- 
strated experimentally in j^J. 

For fixed cavity parameters solitons can exist either for 
A continuously varying in a certain range or for only a 
discrete set of A values. We will show that the latter situ- 
ation is realised for 7 S) j ^ 0. We thus assert that the cav- 
ity selects the frequency difference between the signal and 
idler when a parametric soliton is excited. Note that the 
related problem of frequency selection has been studied 
in the general context of the complex Ginzburg-Landau 
equation (CGLE) with subcritical bifurcation jl6) , which 
also can be applied to describe laser with saturable ab- 
sorber [ p^yT^ ]. The possibility of the approximate re- 
duction of Eqs. (1) to the CGLE in the different limits 



has been demonstrated in M. However, the problem of 
the frequency selection by solitons in OPO has not been 
previously formulated even within the framework of the 
CGLE approximation. This fact has in turn prevented 
construction of the family of stationary soliton solutions 
of Eqs. (2) and rigorous study of their stability. 

To start our analysis of the frequency selection, prob- 
lem, we define the energy parameters Q m = J dx\E m \ 2 
and the energy imbalance (Q_ = Q s — Qi). Manipulation 
of Eqs. (1) reveals that the rate of change of Q- is given 
by 

d t Q- = -2 7 +0- - 2 7 _Q+ = 27.Q, - 2 7s Q s , (3) 

here j± = (j s ± 7i)/2 and Q + = Q s + Qi- Thus for any 
steady state solution, such as soliton solutions of Eqs. 

(2) , the condition 

7s Qs = liQi, (4) 

must be satisfied. Eq. ([|) is consistent with the expec- 
tation that the field with the smaller losses will have the 
larger energy. 

To further interpret relations (||), (^) let us recall 
that in free propagation downconversion is a Hamilto- 
nian process. Then by Noether's theorem every con- 
tinuous symmetry implies a corresponding conservation 
law, see e.g. (TTJ . Cavity losses destroy the Hamiltonian 
structure of the problem, see Eqs. (1), and the input 
pump breaks the phase symmetry linked with conserva- 
tion of the total energy (Q s + Qi + 2Q p ). The symmetry 
(E g ,Ei) -> (E s e^ , Eie-W) in the differential phase V, 
however, survives in the cavity configuration, so how do 
the losses change the associated law dtQ- = 0? Self- 
evidently the relation (^) becomes this conservation law 
in the Hamiltonian limit, which suggests a more general 
link between this energy relation and the symmetry in 
the differential phase. 

Now using cavity solitons as an example we will demon- 
strate that condition (||) constrains the frequency differ- 
ence of the signal and idler components of the intracavity 
field. We consider signal and idler losses small compare 
to the cavity detunings 7 s ,i/|<5| <C 1 and set 6=1. Con- 
servation of Q- is restored for -f St i = and thus Eq. (^) 
is satisfied automatically. A soliton family then exists 
for 5 < 0, with A continuously varying in the interval 

—5). Outside this interval exponential localization of 
A s ^ is not possible. Now suppose that for 7 S ~ e -C 1 
A becomes a slow function of t, i.e. <9(A ~ e, then Eq. 

(3) immediately gives an equation for the adiabatic evo- 
lution of A: d t A ■ d A Q- = -7+Q_(A) - 7 _Q + (A). 
For stationary soliton solutions dtA — and intersec- 
tions of the curve Q s /Qi as a function of A with the line 
Qs/Qi = li/ls give selected values of A. We plot in Fig. 
1(a) the existence curve corresponding to the numerically 
built Jig] single-hump soliton family in the plane (/z, A) 
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for S = — 1 and losses of the signal and idler of order sev- 
eral percent. Dots in Fig. 1(a) correspond to values of 
A obtained by the perturbative method. The agreement 
is excellent, which also indicates that the limit of small 
cavity losses is non-singular. The latter is a necessary 
condition for a linkage between Eq. (3) and the differ- 
ential phase symmetry. Typical transverse profiles of the 
soliton components are presented in Fig. 1(b). 

To study stability of the solitons we seek solu- 
tions of Eqs. (|l|) in the form (A m (x) + e(u m (x,t) + 
iw m (x, t)))e inmt . After standard linearization we derive 
<9t£ = ££, where £ = (u s ,Ui,Up,w s ,Wi,w p ) T and C is 
the linear non-sclf-adjoint differential operator. The dis- 
crete spectrum of C has been found numerically using 
second-order finite differences. Any discrete eigenvalue 
of C with positive real part makes the soliton unstable. 
The new soliton family turns out to be stable over the 
section {A, H) in Fig. 1(a). The Hopf instability of the 
(H, B) branch, and the instability of the (O, A) branch 
are similar to the case of the degenerate OPO |l(|. At 
the point B the single-hump branch bifurcates back in ji, 
initiating a sequence of multi-hump solitons. 

The above perturbation approach to find selected val- 
ues of A requires 7 Sj i small, which is satisfied in most 
practical situations, but we also assumed 8 ~ O(l). This 
fails if 7s,i/|5| ~ 0(1) or ~ e" 1 , which physically means 
that the cavity becomes tuned close to resonance with 
the signal and idler fields. Then terms proportional to 5 
in Eqs. (2) should also be considered as perturbations. 
In this case Eqs. (2) simply do not have solitary solutions 
in zero order and therefore Q St i can't be considered func- 
tions of A. To overcome this difficulty we used another 
perturbation approach, also based on the ^-symmetry. 
Note first that the cavity solitons become wider on ap- 
proach to resonance and so to avoid large compu- 
tational windows it is convenient to take large &, e.g. 
b = 2/T s . Then 7 m ~ 1 physically corresponds to small 
losses. Now, observing that Ai — A s is a solution of Eqs. 
(2) for 7 S = 7, and A = 0, we assume 7_, |A| ~ e. At 
first order in e 



(5) 



where = (— ImA s , ImA s , 0, ReA s , — ReA s , 0) T is the 
neutral mode generated by the ^-symmetry, i.e. = 0, 
and V = (ReA s ,-ReA S7 0,ImA s ,~ImA S7 0) T . Eq. (5) 
immediately yields 



A = 7- 



(C^ICv) 



(6) 



where the new vector £^ is the neutral mode of C\ i.e. 
C^C-ip = 0, generated by the ^-symmetry (which can be 
found numerically) . As usual, (..|..) defines inner product 
in L,2- We again find an excellent agreement between Eq. 
(0) and numerical solutions of Eqs. (2), see Fig. 2(a). 



In this case, stability analysis reveals that the solitons 
are stable along the entire segment (A, B), with B again 
marking the transition to multi-hump soliton solutions. 
We also found that CP |Cv> )/(&/>!&/>) is positive throughout 
a wide region of parameters. This implies that the sign of 
7_ determines the sign of A, at least for parameter val- 
ues where our perturbative approach is valid. Fig. 1(a) 
and selective numerical checks in parameter regions far 
outwith the validity of our perturbative methods sup- 
port a conjecture that the relationship sgn^^ = sgnA 
has a general character. Note finally that any difference 
in diffraction constants a s and a, will also affect the fre- 
quency selection, but we leave this mechanism for future 
analysis. 

To assess possibilities of experimental observation of 
these cavity solitons let us take as an example a 1cm long 
monolithic planar waveguide cavity with — 20pm/ V, 
typical for a noncritically phase matched potassium nio- 
bate crystal. Suppose the waveguide to be ~ 1mm wide 
and ~ 1/xm thick, excited by an elliptical pump beam at 
frequency ~ l0 15 Hz. Other parameters as for Fig. 2(a) 
then imply the following real world values: pump power 
~ /i 2 x 1W, selected frequency difference ~ il x 10 9 Hz 
and cavity soliton width ~ 30/xm. 

To excite cavity solitons one can apply spatially lo- 
calised optical pulses at any of the three frequencies. Op- 
timisation of the pulse parameters goes beyond of our 
present scope. However, to demostrate how the selection 
of the frequency difference happens, we show in Fig. 3 
results of direct simulation of Eqs. (0) with, as initial 
condition, a pulse at frequency oj s with duration around 
0.1t c , intensity several times the peak soliton intensity, 
and width of order the soliton width. Intensities of all 
three components of the excited soliton become constant 
after a transient period, see Fig. 3(a), while the real 
parts of the signal and idler fields exhibit undamped os- 
cillations, confirming selection of A with the predicted 
value, see Fig. 3(b). 

The main physical conclusion which can be drawn from 
the above results is that, if the signal has the smaller 
losses its frequency is pulled away from the cavity res- 
onance while the idler frequency is pulled towards reso- 
nance and vice versa, see Figs. 1(a), 2(a). If the signal 
and idler losses are equal then the selected value of A 
is zero, see Eq. (||), which implies that both the idler 
and the signal are equally detuned from cavity resonance. 
Thus the cavity structure balances the energies of the 
idler and signal components during soliton excitation, in 
accord with Eqs. (||), (Q). The understanding of this 
frequency selection mechanism has allowed us to recon- 
struct an entire family of single-hump cavity solitons and 
to study their stability. Our results are likely to find ap- 
plications in interpretation of other spatio-temporal phe- 
nomena in nondegenerate OPOs and also to be relevant 
in other intracavity parametric processes with symmetry 
in the differential phase, e.g., second harmonic generation 
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with competing parametric process p0[ and nondegener- 
ate four- wave mixing |21|. 
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FIG. 1. (a) Existence curve in the (fi, A)-plane for the fam- 
ily of single-hump cavity solitons, demonstrating selection of 
A. Full line corresponds to the results of numerical contin- 
uation, dotted line was obtained using perturbative methods 
(see text) in the limit 7 S ,;/|<5| < 1: 8 — 8 P = — 1, j s = 0.04, 
ji = 0.05, 7 P = 0.1, b — 1. Solitons are stable only in the 
small interval (A,H), see text for details, (b) Soliton trans- 
verse profile for \i = 0.3, A = -0.1255007. Full line: \A S \, 
dashed line: \Ai\, dotted line: \A P \. 
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FIG. 2. (a) Existence curve in the (/x, A)-plane for family 
of single-hump cavity solitons, demonstrating selection of A. 
Full line corresponds to the results of numerical continuation, 
dotted line was obtained using Eq. (6): 8 — —1.8, 5 P = —4, 
7 S = 1, 7, = 0.95, 7 P = 2, b = 2/T s . Solitons are stable only 
in the interval (A,B). (b) Plots showing dependence of Q m 
vs fi. Parameters as for (a). Full line: Q 3 , dashed line: Qi, 
dotted line: Q p . 
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FIG. 3. Soliton excitation by a localized signal-field pulse. 
Temporal evolution of (a) \E m \ and (b) ReE m at x = 0. Full 




line: signal, dashed line: idler, dotted line: pump. Other pa- 
rameters as for Fig. 1(b). Our predicted A corresponds to a 
period ~ 50 time units. 
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